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Abstract 

We successfully synthesized the zinc-verdazyl complex [Zn(hfac) 2 ]-(o-Py-V) [hfac = 1,1,1,5,5,5- 
hexafluoroacetylacetonate; o-Py-V = 3-(2-pyridyl)-l,5-diphenylverdazyl], which is an ideal model 
compound with an S = 1/2 ferromagnetic-antiferromagnetic alternating Heisenberg chain (F-AF 
AHC). Ab initio molecular orbital (MO) calculations indicate that two dominant interactions Jp 
and Jaf form the S = 1/2 F-AF AHC in this compound. The magnetic susceptibility and magnetic 
specific heat of the compound exhibit thermally activated behavior below approximately 1 K. 
Furthermore, its magnetization curve is observed up to the saturation field and directly indicates 
a zero-field excitation gap of 0.5 T. These experimental results provide evidence for the existence 
of a Haldane gap. We successfully explain the results in terms of the S = 1/2 F-AF AHC through 
quantum Monte Carlo calculations with |Jaf/Jf| = 0.22. The ab initio MO calculations also 
indicate a weak AF interchain interaction J' and that the coupled F-AF AHCs form a honeycomb 
lattice. The J' dependence of the Haldane gap is calculated, and the actual value of J' is determined 
to be less than 0.01 \Jp\- 

PACS numbers: 75.10.Jm 
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I. INTRODUCTION 


Quantum spin systems exhibit unique many-body phenomena caused by strong quan¬ 


tum fluctuations. According to Haldane’s conjecture of 1983 [l|, Heisenberg antiferromag¬ 
netic (AF) chain with integer spin values has a finite energy gap between the non-magnetic 
ground state and the first excited state (Haldane gap), whereas that with half-integer spin 
values has no energy gap. This conjecture has stimulated investigations on Heisenberg 
AF chains, which established the presence of a Haldane gap for integer spin values Oft 
The Haldane state in a Heisenberg AF chain with integer spin values is described as a 
valence-bond-solid (VBS) j^, 7], where each integer spin is considered as two half-integer 
spins forming a singlet state between different sites, and the AF correlation function decays 
exponentially. The Haldane state is characterized by the string order parameter defined 
by ° z str = lim | i—j\—¥oo ( -<Si exp (in Em Sf)Sj>) This parameter indicates hidden 

topological long-range order (LRO) and has a specific value in the Haldane state {lO, 111 ]. 

The valence-bond picture of the Haldane state can be mapped onto a strong ferromagnetic 
(F) coupling limit of an F-AF alternating Heisenberg chain (AHC) with half-integer spin 
values. Therefore, the 5 — 1/2 F-AF AHC has been intensively studied in relation with the 
Haldane state in the 5 = 1 Heisenberg AF chain. The thermodynamic properties of the 5 = 
1/2 F-AF AHC were first calculated by Borras-Almenar et al. nearly two decades ago 12]. 
The field-induced phase transitions to a Luttinger liquid (LL) belong to the same universality 
class as the 5=1 Heisenberg AF chain for both 5 = 1/2 F-AF and AF-AF AHCs, and the 


former and t 
respectively 


re latter have attractive and repulsive interactions between spinless fermions, 


13]. The dynamical structure factor indicates that the first excited state is 


located at q = 7t/2 for the 5=1/2 F-AF AHC similarly to that for the 5 = 1 Heisenberg 
AF chain, while that for 5 = 1/2 AF-AF AHC is located at q = 0 fl4. 15]. Furthermore, the 
neutron scattering intensity is predicted to show a peak at an incommensurate position for 
the 5 = 1/2 F-AF AHC 14]. Hida investigated the ground-state properties of the 5 = 1/2 F- 
AF AHC for various ratios of exchange constants 16| Il8( . There is no discontinuous change 
in the ground state associated with a phase transition between the Haldane (| ,/p| Jaf) 
and AF dimer (| Jp| -C Jaf) phases. The energy gap and string order parameter have finite 
values for all the ratios of exchange constants and change continuously from the Haldane 
phase to the AF dimer phase. It is suggested that the string order parameter represents the 
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strength of the localized singlet correlation, which must be useful to distinguish the VBS 
from other disordered states. 

From the experimental point of view, much effort has been directed at realizing 5 — 1/2 
F-AF AHCs, and several candidates have been reported thus far. Some of the investigated 


compounds indicate the disappearance of the Hale 
due to relatively strong interchain interactions 19 


ane gap and a phase transition to LRO 


20 |. The copper-based IPA-CuCb was 


considered the ideal example until inelastic neutron scattering (INS) experiments were per¬ 


formed 


22], In fact, the magnetic system in IPA-CUCI 3 is characterized as a spin ladder 


with strongly coupled ferromagnetic rungs 


22 ], Although an exited energy gap is observed 


in DMA-CUCI 3 through INS experiments, the effective LRO of spins, excluding those form¬ 


ing F-AF AHCs, renders the intrinsic Haldane-like behavior indistinct 


23 


251] . At present, 


a-CuNb 206 and NasCuoSbOe are promising candidates with 5=1/2 F-AF AHCs because 
their magnetic susceptibility, NMR, and INS measurements demonstrate Haldane-like be¬ 


havior 


26 


30]. However, the entire phase diagram extending to the gapless phases of these 


materials have not been experimentally reported, most likely owing to large energy gaps 
requiring high magnetic fields. The entire range of the magnetization curve up to the sat¬ 
uration field, which describes not only the Haldane gap but also the filling of a fermionic 
band in the field-induced LL phase, is essential to further our understanding of the quantum 
effect in the crossover from Haldane to AF dimer phases. 

In this paper, we successfully synthesized the zinc-verdazyl complex [Zn(hfac) 2 ]-(o-Py-V) 
[hfac = 1,1,1,5,5,5-hexafluoroacetylacetonate, o-Py-V = 3-(2-pyridyl)-l,5-diphenylverdazyl], 
which is an ideal model compound with the 5 = 1/2 F-AF AHC, and investigated its 
crystal structures at room and low temperatures. Ab initio molecular orbital (MO) calcula¬ 
tions indicated that two dominant interactions form an 5 = 1/2 F-AF AHC. We observed 
thermally activated behavior associated with an energy gap in the magnetic susceptibility 
and the magnetic specific heat. The magnetization curve is observed up to the saturation 
field and directly indicates a zero-field excitation gap of 0.5 T. We successfully explained 
these experimental results in terms of the 5 = 1/2 F-AF AHC through quantum Monte 
Carlo (QMC) calculations. Furthermore, we considered the effects of weak AF interchain 
interactions forming an 5 = 1/2 honeycomb lattice. 
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II. EXPERIMENTAL 


We synthesized o-Py-V through a conventional procedure 31]. A solution of 
[Zn(hfac) 2 ]-2H 2 0 (176 mg, 0.34 mmol) in 20 ml of heptane was refluxed for 10 min at 
70 °C. A solution of o-Py-V (107 mg, 0.34 mmol) in 10 ml of CH 2 C1 2 was slowly added, 
and stirring was continued for 1 h. After the mixed solution cooled to room temperature, 
a dark-green crystalline solid was separated by filtration and washed with pentane. The 
dark-green residue was recrystallized using CH 2 C1 2 in an acetonitrile atmosphere (155 mg, 
0.195 mmol), after which the yield was approximately 57.4%. 

X-ray intensity data were collected using a Rigaku AFC-7R mercury CCD diffractometer 
and a Rigaku AFC-8R mercury CCD RA-micro7 diffractometer at 293 and 25 K, respec¬ 


tively, with graphite-monochromated Mo Ka radiation and Japan Therma. 


XR-HR10K. The structure was solved by a direct method using SIR2004 32] and was re 


fined with SHELXL97 


Engineering 


33]. The structural refinement was carried out using anisotropic and 


isotropic thermal parameters for the nonhydrogen atoms and the hydrogen atoms, respec¬ 
tively. All the hydrogen atoms were placed at the calculated ideal positions. The magnetic 
susceptibility and magnetization curves were measured using a commercial SQUID mag¬ 
netometer (MPMS-XL, Quantum Design) and a capacitive Faraday magnetometer with a 
dilution refrigerator. The experimental results were corrected for the diamagnetic contri¬ 
bution of -3.29xlCD 4 emu-moh 1 , which was calculated using Pascal’s method. The specific 
heat was measured with a commercial calorimeter (PPMS, Quantum Design) by using a 
thermal relaxation method above 2.0 K and by using an adiabatic method between 0.35 K 
and 2.0 K. All experiments were performed using small randomly oriented single crystals 
with typical dimensions of 2.Ox 1.0x0.5 mm 3 . 

Ab initio MO calculations were performed using the UB3LYP method as broken- 
symmetry (BS) hybrid density functional theory calculations. All calculations were per¬ 
formed using the Gaussian 09 program package and 6-31G basis sets. The convergence 
criterion was set at 10 -8 hartree. For the estimation of intermolecular magnetic interac¬ 
tions, we applied an evaluation scheme to previously studied multispin systems using the 
Ising approximation 34]. 

The QMC code is based on the directed loop algorithm in the stochastic series expansion 
representation 35]. The calculations for the S — 1/2 F-AF AHC was performed for N = 256 
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under the periodic boundary condition, where N denotes the system size. The QMC code for 
the evaluation of spin gap is based on the multi-cluster loop algorithm in the continuous-time 
path-integral representation. The spin gap is calculated by the second-moment method; that 
is, the gap is the inversed correlation length in the imaginary-time direction. All calculations 
were carried out using the ALPS application 2.1361. |37|. 


III. RESULTS 


A. Crystal structure and magnetic model 


The crystallographic data are summarized in Table I 


38], and the molecular structure is 


shown in Fig. 1(a). The verdazyl ring (which includes four nitrogen atoms), the two upper 
phenyl rings, and the bottom pyridine ring are labeled as Ri, R 2 , R3, and R 4 , respectively. 
The MO calculation indicates that approximately 58% of the total spin density is present on 
Ri. While R 2 and R 3 each account for approximately 17 % of the relatively large total spin 
density, R 4 accounts for less than 6 % of the total spin density. Since Zn(hfac ) 2 has a low spin 
density, it works as a spacer between verdazyl radicals, resulting in the low dimensionality of 
the magnetic lattice. We focus on the structural features related to the o-Py-V to consider 
intermolccular interactions. We evaluated the intermolecular magnetic interactions of all 
molecular pairs within 4.0 A at both 293 and 25 K through the ab initio MO calculations. 
Consequently, we found that there are three types of dominate interactions related to M 0 - 
Mi, Mo-M 2 , and M0-M3 molecular pairs, as shown in Figs. 1(b)-(d). The Mo-Mi molecular 
pair has an N-C short contact di, which is doubled by inversion symmetry, as shown in Fig. 
1 (b), and it is approximately 3.48 A at both 293 and 25 K. A strong F interaction Jp is 
evaluated between these molecules. The M 0 -M 2 molecular pair has a C-C short contact d 2 , 
which is also doubled by inversion symmetry, as shown in Fig. 1(c). The values of d 2 at 
293 and 25 K are 3.48 and 3.39 A, respectively. A relatively strong AF interaction Jaf is 
evaluated between these molecules. The Mo-Mi and Mo-M 2 molecular pairs are alternately 
aligned along the c-axis, as shown in Fig. 1(e), and form an S = 1/2 F-AF AHC consisting 
of the Jp and the J\p, as shown in Fig. 1(f). The M 0 -M 3 molecular pair, which is related by 
inversion symmetry, is associated with interchain interaction J' and has a C-C short contact 
d 3 . The values of d 3 at 293 and 25 K are 3.51 and 3.44 A, respectively. Although there is 
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relatively short contact, the small overlap of the n orbitals, which expand perpendicular to 
the planes, should render J' weak. This interchain interaction forms a honeycomb lattice in 
the ac-plane. The evaluated values of the exchange interactions are summarized in Table II. 
They are defined with the Heisenberg spin Hamiltonian: 


H = S 2i+i + Jaf + J'Y.Sk-Su 


( 1 ) 


i i kl 

where S is an 5=1/2 spin operator. It is remarkable that the absolute value of the interchain 
interaction J' at 25 K has a small value, while the other interactions are more dominant 
at 25 K than at 293 K. Thus, the one-dimensionality is expected to be enhanced in the 
low-temperature regions, which should provide an intrinsic magnetic behavior reflecting the 
ground state of the 5 = 1/2 F-AF AHC. 


B. Magnetic and thermodynamic properties 


Figure 2 shows the temperature dependence of the magnetic susceptibility (y = M/H) 
at 0.1 T (H = 1000 Oe). We observed a broad peak at approximately 1.0 K, below which y 
decreases with decreasing temperatures, indicating the existence of a nonmagnetic ground 
state separated from the excited states by an energy gap. Above 100 K, it follows the 
Curie-Weiss law, and the Weiss temperature is estimated to be $w = +3.0(3) K, which 
indicates the dominant contribution of Jp. We evaluated the paramagnetic impurities to 


be approximately 1.2% of all spins, which is defined to fit the fol 
and is close to those evaluated in other verdazyl radical crystals 


owing calculated result 


39 


41], and subtracted 


it from the raw data by assuming conventional paramagnetic behavior, as shown in Fig. 
2. The contribution of Jp appears in the temperature dependence of yT, which increases 
with decreasing temperatures down to approximately 12 K, as shown in the inset of Fig. 
2. Below 12 K, yT decreases with decreasing temperature, indicating the existence of AF 
interactions. 

Figure 3 shows the magnetization curve at 0.13 K. We subtracted a small paramagnetic 
contribution given by the Brillouin function in the low-field region that is evaluated to 
account for 1.2% of all spins. The saturation value of 0.98 /Xp/f-u. indicates that the purity 
of the radicals is approximately 98 %, and we consider this purity in the following analysis. 
The corrected magnetization curve without the paramagnetic contribution clearly exhibits 
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TABLE I: Crystallographic data for [Zn(hfac) 2 ]-(o-Py-V) 


Formula 

C29H18F 

i2N 5 04Zn 

Crystal system 

Monoclinic 

Space group 

P2ifc 

Temperature (K) 

293(2) 

25(2) 

Wavelength (A) 

0.7107 

a(A) 

9.130(5) 

8.8144(17) 

6 (A) 

31.839(16) 

31.582(6) 

C (A) 

10.885(6) 

10.826(2) 

(3 (degrees) 

93.286(5) 

92.772(4) 

v (A 3 ) 

3159(3) 

3010.3(10) 

Z 


4 

£>calc (g cm -3 ) 

1.669 

1.752 

Total reflections 

4796 

4390 

Reflection used 

4289 

3991 

Parameters refined 

460 

R [I > 2<t(/)] 

0.0661 

0.0787 

R w [I > 2 a(I)} 

0.1669 

0.1922 

Goodness of fit 

1.065 

1.029 

CCDC 

1020097 

1020098 


TABLE II: Evaluated magnetic interactions through the ab initio MO calculations at 293 and 25 
K and the experimental analysis by using QMC calculations 


Method \Jaf/Jf\ \J'/Jf\ Jf/Izb (K) JafAb (K) J'/k B (K) 


MO at 293 K 

0.29 

0.13 

-17.2 

5.0 

2.3 

MO at 25 K 

0.34 

0.05 

-21.2 

7.2 

1.1 

QMC analysis 

0.22 

< 0.01 

-12.8 

2.8 

< 0.1 











FIG. 1: (color online) (a) Molecular structure of [Zn(hfac) 2 ]-(o-Py-V). Molecular packing of neigh¬ 
boring (b) Mo-Mi, (c) M 0 -M 2 , and (d) M 0 -M 3 molecular pairs, which are associated with Jf, Jaf, 
and J ', respectively. Broken lines indicate C-C and N-C short contacts. Hydrogen atoms are 
omitted for clarity, (e) Crystal structure forming an alternating chain consisting of Jp and Jaf 
along the c-axis, and (f) the corresponding S = 1/2 F-AF AHC. 


a small zero-field excitation gap of approximately 0.5 T, as shown in Fig. 3 and its upper 
inset. The lower inset of Fig. 3 shows the field derivative of the corrected magnetization 
curve ( dM/dB ). The steep increase of dM/dB at the critical held region evidences the 
gapped behavior. This gapped behavior of the magnetization curve strongly suggests the 
realization of a Haldane state in the 5 = 1/2 F-AF AHC. 

The inset of Fig. 4 shows the experimental result of the total specific heat C p at 0 and 
1.0 T. The magnetic specific heats C m are obtained by subtracting the lattice contribution 
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FIG. 2: (color online) Temperature dependence of magnetic susceptibility (y = M/H ) of 
[Zn(hfac) 2 ]-(o-Py-V) at 0.1 T. The open circles denote raw data, and the open triangles are cor¬ 
rected for the paramagnetic term due to the impurity. The inset shows the temperature dependence 
of \T. The solid lines with open triangles, squares, and circles represent the calculated results for 
the S = 1/2 F-AF AHC with a = |Jaf/Tf| = 0.10, 0.22, and 0.40, respectively. 

by assuming Debye’s T 3 -law as 0.020T 3 (J/mol K), which corresponds to about 23 K of the 
Debye temperature. There is no sharp peak associated with a phase transition to an ordered 
state. A clear Schottky-like peak associated with an excitation gap below approximately 1.0 
K at zero field is observed, which is consistent with the existence of a Haldane gap. 


IV. DISCUSSION 

We discuss the ground state of the present spin model. The MO calculations show that 
the dominant intermolecular interactions Jp and Jaf form the S — 1/2 F-AF AHC, and 
the one-dimensionality is enhanced in the low-temperature regions, as summarized in Table 
II. Accordingly, we analyzed the experimental results in terms of the S' = 1/2 F-AF AHC 
as shown in Fig. 1(f). We calculated the magnetic susceptibility, magnetization curve, and 
magnetic specific heat as a function of a = \Jaf/Jf\ by using the QMC method. As is often 
the case with conventional radical compounds, we assume the Heisenberg spin Hamiltonian 
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FIG. 3: (color online) Magnetization curve of [Zn(hfac) 2 ]-(o-Py-V) at 0.13 K. The open triangles 
denote raw data, and the values of the vertical axis have been shifted by 0.1 hb/ f.u. for clarity. 
The closed circles are corrected for the paramagnetic term due to the impurity. The solid line with 
open squares represents the calculated result for the 5 = 1/2 F-AF AHC with a = |Jaf/«/f| = 
0.22. The upper and the lower insets show the expansion of the low-field region and held derivative 
dM/dB , respectively. 


and g = 2.00. Although there is no distinct dependence of x on a, the maximum value of 
X T is very sensitive to the change in a, as shown in the inset of Fig. 2. We obtained good 
agreement between the experiment and calculation by using the parameters Jf/^b = -12.8 K 
and JafAb = 2.8 K (a = 0.22), as shown in Fig. 2 and its inset. The slight difference, espe¬ 
cially in the low-temperature region, might originate from the inaccuracy in the subtraction 
of the paramagnetic impurity. In addition, the calculated result corresponds to the limiting 
behavior H —> 0, which also causes a slight difference reflecting the change in the energy gap 
at the experimental field of 0.1 T. A comparison between the experimental results and ab 
initio calculations on verdazyl-based materials 4lJ, [42] shows that the obtained parameters 
are consistent with those evaluated from the MO calculation. The magnetization curves 
calculated using the QMC method also well reproduced the experimental results using the 
same parameters, as shown Fig. 3 and its insets. Because the energy scale of the exchange 
interactions is relatively small, the finite temperature effect is enhanced in the present sys- 
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FIG. 4: (color online) Temperature dependence of magnetic specific heat of [Zn(hfac) 2 ]-(o-Py-V) 
at 0 and 1.0 T. The solid lines with open circles and triangles represent the calculated results for 
the 5 = 1/2 F-AF AHC with a = |Jaf/«^f| = 0.22 at 0 and 1.0 T, respectively. The inset shows 
the total specific heat. 


tem, resulting in the finite slope in the gapped region below 0.5 T and the gradual change 
at the critical held. We find a slight deviation in the vicinity of the saturation held. In the 
experimental result, the nonlinear behavior originating from the one-dimensional (ID) quan¬ 
tum fluctuation is marginally suppressed by the weak interchain interactions forming the 
two-dimensional (2D) honeycomb lattice, resulting in the difference between the experiment 
and calculation. The calculated result of the magnetic specihc heat also well reproduces 
the Schottky-like peak at 0 T, as shown in Fig. 4. Considering the gapless phase above 
0.5 T, the broad peak observed at 1.0 T indicates the development of ID short-range order. 
These different origins of the broad peaks indeed give rise to a large difference in the peak 
shapes and values. The difference between the experiment and calculation should arise from 
the inaccuracy in the subtraction of the lattice contribution, which we assumed largely as 
Debye’s T'Maw. The ground state for the S = 1/2 F-AF AHC with a = 0.22 is a Haldane 
state with the excitation gap A ps 0.67 K to the lowest triplet states, which is consistent 
with a previous theoretical study 16]. The corresponding string order parameter is pre- 
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dieted to show an intermediate value between the S — 1 Haldane and S — 1/2 AF dimer 
phases jl6] . The gapless phase is described as an LL and expected to have attractive spinless 
fermions 13]. The present small excitation gap and easily accessible saturation held will 
enable quantitative tests of such properties through various measurements. 

Finally, we consider the contribution of J' quantitatively. The slight deviations in the 
low-temperature regions of the magnetic susceptibility and in the vicinity of the zero and 
the saturation fields of the magnetization curve may indicate slight contributions from the 
weak interchain interaction J'. The S — 1/2 F-AF AHCs coupled by the AF J' form an 
S — 1/2 honeycomb lattice, as shown in Fig. 5. We calculated the value of the energy 
gap A as a function of [3 — \J'/J F \, assuming a = 0.22 and J F /k B = -12.8 K using the 
QMC method at \T/ Jp\ = 0.001. The inset of Fig. 5 shows the system size dependence of 
A for each value of f3. Assuming an energy gap in a 2D system, finite-size data are fitted 
using Atv = Aqo + C'exp(—DA 1//2 ), where C and D are fitting parameters 43 


44J. It is 


clear that the fitting curves are appropriate for j3 < 0.02, as shown in the inset of Fig. 
5. For [3 > 0.03, on the other hand, it is difficult to estimate the value of the energy gap 
in the thermodynamic limit by using the above fitting equation because there is no convex 
downward behavior in the fini te-size data up to N— 576. Therefore, we roughly evaluated the 
energy gap in the thermodynamic limit using a linear fit with the values of the four largest 
system sizes (N =324, 400, 484, and 576) for f3> 0.03. The energy gap decreases dramatically 
with increasing (3 and almost disappears at (3 — 0.04, resulting in the phase transition to an 
AF ordered state, as shown in Fig. 5. The experimental result of the magnetization curve 
shows that A fa 0.5 T (fa 0.67 K); the experimental and calculated results are consistent. 
Therefore, even if the interchain interaction is considered, the value of A should be more 
than 0.4 T (fa 0.54 K), as deduced from the gapped behavior of the magnetization curve. 
Considering that the calculated value of A for (3 = 0.01 is approximately 0.51 K, the actual 
value of f3 is roughly estimated to be less than 0.01. 


V. SUMMARY 

We succeeded in synthesizing a zinc-verdazyl complex [Zn(hfac) 2 ]-(o-Py-V). The ab initio 
MO calculations indicate the existence of dominant interactions J F and Jaf forming an 
S' = 1/2 F-AF AHC in this compound. Thermally activated behaviors are observed in 
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FIG. 5: (color online) /3 = \J'/Jp\ dependence of the energy gap A for the weakly coupled 5 = 1/2 
F-AF AHCs with a = | Jay/Jf\ = 0.22 and Jp/kp = -12.8 K. The inset shows the 1/N dependence 
of |A/Jf| for each /3. The broken lines are fitting curves described in the text. The illustration 
shows the 5 = 1/2 honeycomb lattice consisting of Jp, Jap, and J'. 


the magnetic susceptibility and the magnetic specific heat, and the magnetization curve 
indicates the existence of an energy gap of 0.67 K. These features are described as the 
magnetic behavior in an S = 1/2 F-AF AHC with the parameters Jp/ks = -12.8 K and 
JAp/kp = 2.8 K (|Jaf/<7f| = 0.22) by using the QMC method. Furthermore, the MO 
calculation also indicated the existence of the interchain interaction J' and that the coupled 
F-AF AHCs form a honeycomb lattice. We consider the contribution of J' to the value of 
the Haldane gap by using the QMC method, and the actual value of |J'/J F | is evaluated 
to be less than 0.01. These results demonstrate that [Zn(hfac) 2 ]-(o-Py-V) is an ideal model 
compound with an S = 1/2 F-AF AHC. Further investigations of dynamical properties will 
yield quantitative information about the intermediate phase between 5 = 1 Haldane and 5 
= 1/2 AF dimer phases and the field-induced LL phase of attractive spinless fermions. 
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